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Abstract 

The ATLAS experiment measures yields of isolated photons and of Z and W bosons via leptonic decay modes in 
Pb+Pb collisions at -^5^=2.76 TeV. The data samples used in the analysis were obtained in the year 2010 and year 
2011 LHC runs and correspond to 5 /jb -1 and 0.15 nb~' of integrated luminosity respectively. The measured yields 
of all bosons are consistent with the scaling proportional to the number of nucleon-nucleon collisions. The transverse 
momentum distributions are measured for isolated photons and for Z bosons for different centrality bins. The shapes 
of measured distributions do not change with centrality. The transverse momentum and rapidity distributions of Z 
bosons are consistent in shape with Pythia model. 
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1 1. Introduction 

2 The study of the heavy ion (HI) collisions carried out by the experiments at the Relativistic Heavy Ion Collider 

3 (RHIC) at BNL, and the Large Hadron Collider (LHC) at CERN, demonstrated that the hot and dense matter produced 

4 in the interactions of heavy nuclei at high energy imposes a significant energy loss on the energetic color charge car- 

5 riers penetrating such a medium [ 1 ]. An understanding of this phenomenon requires measuring the initial production 
s rates of the particles in HI collisions before they lose energy. The best candidates to perform such measurements 

7 are particles which do not interact via the strong force, such as W, Z bosons and photons. The PHENIX experi- 

8 ment at RHIC measured the properties of highly energetic photons [2]. The ATLAS and CMS experiments at LHC 

9 provided the first measurements of Z, W and photons at the LHC energy [3] 01 [5] |6] |7) . Due to the complexity of 

10 these measurements and limited statistics of the first LHC Pb+Pb run, the results have relatively large uncertainties. 

11 Within their uncertainties the measurements show that in HI collisions the production rates of particles at high mo- 

12 mentum transfer are proportional to average nuclei thickness function ((7aa)X commonly expressed also as mean 

13 number of binary collisions, (N co \i), between the nucleons of the colliding nuclei. These two parameters are related 
u as (N co n) = (Taa) X o~p P , where o- pp is the total inelastic cross section of p + p interactions. The calculations of these 
15 parameters are carried out within the framework of the Glauber model 00 IS]. 

is This proceeding presents the latest results from the ATLAS experiment at LHC on the measurement of W, Z 

17 bosons and isolated photon production in Pb+Pb collisions at Vsaw=2.76 TeV. The results measured in different 

18 centrality bins are compared to each other and also to predictions. 



/ Nuclear Physics A 00 (2012) 1^8\ 



2 



2. ATLAS detector 

The ATLAS detector 1 10] at the LHC covers nearly the entire solid angle around the collision point. It consists of 
an inner tracking detector surrounded by a thin superconducting solenoid, electromagnetic and hadronic calorimeters, 
and a muon spectrometer incorporating three large superconducting toroid magnet systems. 

The inner-detector system (ID) is immersed in a 2 T axial magnetic field and provides charged particle tracking in 
the pseudorapidity range \r]\ < 2.f[j] The high-granularity silicon pixel (Pixel) detector covers the vertex region with 
three layers of pixels of 400 x 50 fim 2 and it is followed by four layers of double sided silicon microstrip tracker (SCT) 
and the transition radiation tracker. 

The calorimeter system covers the range \t]\ < 4.9. Within the region \rj\ < 3.2, electromagnetic calorimetry is 
provided by barrel and end-cap high-granularity lead-liquid argon (LAr) calorimeters, with an additional thin LAr 
presampler covering \rj\ < 1.8. 

The muon spectrometer (MS) comprises separate trigger and high-precision tracking chambers measuring the de- 
flection of muons in a magnetic field generated by superconducting air-core toroids. The precision chamber system 
covers the region \rj\ < 2.7 with three layers of monitored drift tubes (MDT), complemented by cathode strip cham- 
bers (CSC) in the innermost layer of the forward region, where the background is highest. The muon trigger system 
covers the range I77I < 2.4 with resistive plate chambers in the barrel, and thin gap chambers in the end-cap regions. 

The two minimum-bias trigger scintillator (MBTS) counters, covering 2.1 < I77I < 3.9 on each side of the nominal 
interaction point register charged particle hits with 16 scintillating pads on each side. Two zero-degree calorimeters 
(ZDC), each positioned at 140 m from the collision point, detect neutrons and photons with \rj\ > 8.3. 



The analysis of W bosons uses 2010 LHC Pb+Pb collision data corresponding to an integrated luminosity of 
approximately 5/ib~'. The analysis of isolated photons and Z bosons is based on the 2011 LHC Pb+Pb collision 
data and corresponds to an integrated luminosity of approximately 0.15nb~'. Both data samples were obtained at the 
energy of ysMV=2.76 TeV. 



3.1. Event selection 

All selected events are required to satisfy the Minimum Bias (MB) event selection. All events in 2010 run were 
taken with the MB trigger, requiring the coincidence of the ZDC signals on both sides, and the event vertex recon- 
structed offline. With increasing event rates in 2011 run, this condition was modified. On a trigger level the MB 
required transverse energy (Ej) deposition of E T >50 GeV in the ATLAS calorimeter or a coincidence of the ZDC 
signals on both sides and a track in the ID system. In offline analysis all events were required to have an event vertex 
reconstructed by the ID tracking system and also an MBTS timing signal coincidence | Af| < 3 ns. The total number 
of sampled events is different in different analyses. It varies from approximately 35 x 10 6 MB events selected for the 
analysis of W in 2010 run to 1.03 x 10 9 MB sampled events in the analysis of the Z — > ee decay channel ifTTl in 201 1 
data. 

3.2. Event centrality 

In HI collisions, "centrality" reflects the overlap volume of the two colliding nuclei, controlled by the impact 
parameter. That overlap volume is closely related to the average number of participants, (N pmt ), the nucleons which 
scatter inelastically in each nuclear collision and the number of binary collisions, (N co n), between the nucleons of the 
colliding nuclei. 

The Pb+Pb collision centrality is measured using the sum of transverse energy (2 £t) deposited in the Forward 
Calorimeter (FCal) over the pseudorapidity range 3.1 < \rj\ < 4.9 calibrated at the electromagnetic energy scale 0~2). 
The fraction of events with more than one collision in the 2010 run was negligibly small, and in the 201 1 run it was 
estimated not to exceed 0.05% independent of centrality except for the very central collisions. A cut on the FCal 
energy of £ £t<3.8 TeV was applied to prevent pileup contamination in the analyzed data. 

1 The ATLAS reference system is a Cartesian right-handed coordinate system, with the nominal collision point at the origin. The anticlockwise 
beam direction defines the positive z-axis, while the positive .r-axis is defined as pointing from the collision point to the center of the LHC ring and 
the positive y-axis points upwards. The azimuthal angle <f> is measured around the beam axis, and the polar angle 8 is measured with respect to the 
z-axis. Pseudorapidity is defined as rj = - ln(tan(6/2)). 



3. Analysis 
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3.3. Reconstruction of the W bosons 

Measurement of the W bosons is based on reconstruction of inclusive muon spectrum. Measurements of the muon 
trajectories from both the ID and MS are combined, resulting in a relative momentum resolution ranging from about 
2% at low momentum up to about 3% at transverse momentum (pj) about 50GeV. Muon tracks are required to have 
at least two hits in the Pixel detector, one of them in the first layer, and six or more hits in the SCT. To improve on 
the match of the muon spectrometer track and the inner detector the momentum measured in the muon spectrometer 
must be within 50% of the corresponding measurement in the ID. Left panel of Fig. [T] shows the uncorrected muon 
spectrum obtained in all events in the 2010 run IfTZll . 




Figure 1: Left: Uncorrected inclusive muon px spectrum from Ref. |5 |. The spectrum is fitted (solid line) with two components: signal W — > pv 
(shaded area) simulated with Pythia in p + p collisions, and a background parametrization (dashed line). The invariant mass distributions of 
Z — > ee (middle) and Z — > pp (right) in data and MC integrated over momentum, rapidity, and centrality from Ref. [11]. The MC is reweighted 
to match the centrality distribution in data and normalized in the region 66 GeV < m;; < 102 GeV (I = e,p) indicated by dashed lines. Numbers of 
counts given in the plot correspond to the same mass region. 

The spectrum shows a steep power law fall up to pj of 30 GeV where the presence of muons from the W decay 
appears prominently. The W production yields are obtained performing fits to the muon transverse momentum spectra 
using two input shapes. The first shape, shown as shaded area in the left panel of Fig. [TJ describes the muon pj 
distribution from W ± — > p ± v, and a second describes the background of exponential shape, obtained from studies of 
cc — > p. + X and bb — > p + X decays in p + p simulations and described in [5|. 

For the present analysis, only muons with pi > 7 GeV are considered, and also muon pairs forming the invariant 
mass > 66 GeV are removed from the analysis of the W decays to reduce the background coming from Z. After 
applying cuts, the total number of W ± candidates extracted from the event sample by fitting the templates is 399+™. 

3.4. Reconstruction of the isolated photons 

Photon candidates were identified at the first trigger level (L 1 ) by a cluster formed by ( A(p x A77) = 0.1x0.1 trigger 
towers of the electromagnetic part of the calorimeter, covering a pseudorapidity range \rj\ < 2.5. A cluster transverse 
energy estimated at LI was required to exceed £r=16GeV. The trigger becomes fully efficient to the photons with 
energy exceeding 20 GeV. 

In order to reconstruct photons and electrons in the context of a heavy ion collision, the large background from 
the underlying event (UE) is subtracted from each event. This is performed in the same way as in the heavy ion jet 
reconstruction, explained in detail in Ref. IfTJI . The UE subtraction takes into account the average energy density in 
the calorimeter modulated by the flow effects in the regions excluding jet candidates. 

As described in Ref. [ 14 1, after applying the shower shape cuts, but before applying isolation cuts there are 6435 
photon candidates with pj > 45 GeV and within \rj\ < 1.3 the (0-80)% centrality sample. Additional kinematic 
constraints on photon candidates are used to restrict the analysis to the region with less material in front of the 
calorimeter and to improve the performance of the photon identification. 
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92 The isolation energy is the sum of transverse energies in calorimeter cells within a cone or a radius R lso — 

93 -\jAr/ 2 + Acf> 2 around the photon direction, excluding the cells associated to the photon itself. Figure [2] shows the 
distributions of Ej for R\ so = 0.3 as a function of collision centrality. The data and MC distributions shown in the 




-20 20 40 60 -20 20 40 60 

E T (R. ao =0.3) [GeV] E T (R. so =0.3) [GeV] 



Figure 2: Distributions of photon isolation energy in a Ri so = 0.3 cone for the three centrality bins in data (black points) and for MC (yellow 
histogram), normalized for negative Ej values from Ref. 1 14 1. The differences at large values of Ej can be attributed to the presence of jet 
contamination in the data, which is not present in the MC sample. 

94 

95 figure grow noticeably wider with increasing centrality, as the UE subtraction is only able to subtract the mean energy 

96 in an interval. The effect is reproduced by the MC, except for a contribution from di-jets which are not present in 

97 simulation. An isolation cut of 6 GeV was used in the analysis. The absolute energy scale of the photons derived from 

98 the Pythia events embedded into Hijing UE is withing 1.5% of the generated value. The energy resolution is better 

99 than 3% in the most central collisions. 

100 The photon yields are extracted using the double sideband technique which divides photon candidates according 

101 to their quality and isolation energy as described in lfl4l . This technique provides an accurate estimation of residual 

102 background in the isolated photon sample and signal leakage in the background region. 

103 3.5. Reconstruction of the Z bosons 

104 The Z bosons are reconstructed via pp and ee decay channels. A data sample for this study was collected by the 

105 ATLAS three-level trigger system [ 15] which selected events with high-px muon and electron candidates in 201 1 data 

106 sample. 

107 High-/?x muon candidates were selected using all three levels of the ATLAS trigger system. The LI muon trigger 

108 searches for patterns of hits consistent with muons at certain pj in the trigger chambers within \t]\ < 2.4. If a muon 

109 p T estimate exceeds 4 GeV, the event is accepted for further processing at the High Level Trigger (HLT). The LI 

110 muon algorithm also identifies Regions of Interest (Rol) within the detector to be investigated by the HLT. At the 

111 HLT, each muon candidate's track parameters are recalculated by including the precision data from the MDT or CSC 

112 in the Rol defined by the previous trigger level. Muon candidates are reconstructed either solely from the MS or using 

113 combined data from the MS and ID systems. In addition to the events selected using the Rol-based muon trigger, an 

114 event filtering trigger was used to identify muons with pj > 10 GeV. The filtering was done for all events satisfying a 

115 LI requirement of the ZDC firing on both sides or an Ej > 10 GeV deposition in the calorimeter. 

no In the Z — > pp analysis, single muons may be reconstructed with varying levels of qualitv lfTBI . High quality muons 

117 are reconstructed in the MS and ID subsystems and have a consistent angular measurement in both subdetectors as 

us well as a good match to the primary vertex. At least one muon in a pair, matched to the trigger, is required to be 

119 of such quality. If the second muon in the pair is fully reconstructed as a track by both subsystems, the minimum 

120 p T threshold is set to 10 GeV on both muons. If the second muon fails this condition, both muons are required to 

121 satisfy px>20GeV. An invariant mass window of 66 GeV < m ff < 102 GeV is used to define oppositely charged 

122 muon pairs as Z boson candidates and same-sign charged pairs as a background estimate. In total, 1223 opposite-sign 

123 candidates and 14 same-sign candidates reconstructed in the Z — > pp channel are shown in the right panel of Fig.[T] 

124 Electron candidates are reconstructed at LI requiring Ej to exceed 14 GeV. The UE subtraction is performed for 

125 photon clusters in the same way as for the photon candidates. 
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126 For the Z — > ee analysis, electron candidates are formed using the standard ATLAS reconstruction algorithm, 
requiring the matching of an ID track to a LAr cluster. Electron selection is limited to \rj\ < 2.5 and both electrons 
are required to have £t>20 GeV. The standard electron identification cuts ifTTll used in the p + p environment are not 
suited to the Pb+Pb environment due to the amount of energy from the underlying event deposited in the calorimeter. 
To address this, a different set of cuts was developed to accommodate the modification of the calorimeter variables 
by the presence of the underlying event. The electron identification cuts used were based on the energy balance 
between the track and EM cluster (Ej p) and shower shape variables. Since almost all electrons with £t>20 GeV are 
reconstructed as jets by the heavy-ion anti-A: t R=0.2 jet finding algorithm lfl3l . a further cut was made requiring that 
the Ej of the electron be at least 53% of the jet Ej. 

For the Z — > ee analysis all electrons found in triggered events were paired with each other, requiring that at least 
one electron in the pair was matched to a trigger object. The opposite-sign pairs within an invariant mass satisfying 
66 GeV < m ee < 102 GeV are accepted as signal Z candidates and the same charged-sign pairs in this window 
are taken as an estimate of the combinatorial background. In total, 772 opposite-sign pairs and 42 same-sign pairs 
reconstructed in the Z — > ee decay channel are shown in Fig. [T] 

3.6. Corrections and uncertainties 

The W boson is measured using the MB event sample, for the measurement of Z boson and isolated photons, 



154 



measured with triggered sample, correction for finite efficiency of the trigger system is applied as explained in sec. 3.5 
In Z — > ee channel and in the measurement of isolated photons these corrections are significantly smaller than 1%. 

Reconstruction and identification efficiency of Z bosons was evaluated using 7 x 10 5 Z — > ee events and 5.3 X 10 5 
Z — > pp events with 66 GeV < ntz < 1 16 GeV and \y z \ < 2.5 produced by Pythia event generator [ 18 1 and embedded 
into Pb+Pb events generated by the Hijing |[T9l event generator. The response of the ATLAS detector to the generated 
particles is modeled using GEANT4 [20|. 

The Z boson distributions are corrected for the losses due to detector acceptance and analysis cuts in the bins of 
pi, Tj and centrality. The correction is larger for the Z — > ee channel, as the electron losses are higher. At the very 
edge of the acceptance region the Z — > ee reconstruction efficiency drops to 10% from more than 50% at midrapidity. 
The corresponding numbers for the Z — » pp measurements are 45% and 70% respectively. 

The main source of systematic uncertainties in the measurement of isolated photons are the photon quality cuts 
and photon isolation criteria. The uncertainties reach 20% for each. In the measurement of the Z bosons the main 
contributions in both decay channels are from the uncertainties of efficiency corrections reaching up to 8%. Because 

155 both channels have different sources of systematic uncertainties they can be averaged taking into account the statistical 

156 and systematic uncertainties. Finally, the uncertainties on the (N co \i) reach more than 10% and they play dominant 

157 role in several results. 



158 3.7. Results 

159 The fully corrected p T distribution of isolated photons measured by the ATLAS experiment is shown in the left 

160 panel of Fig. [3] The CMS Collaboration measured production of isolated photons Q in heavy ion collisions at 

161 V?Aw=2.76TeV in the rapidity region \r(\ < 1.44 and an isolation condition of at most 5 GeV transverse energy in a 

162 cone of radius R lso = 0.4. The measurement extends to 80 GeV. The CMS (0-10%) data are superimposed on the most 

163 central bin as well as the CMS p + p data are superimposed on the most peripheral bin. Within the range of both 

164 measurements the results are consistent. 

165 For both Z — > ee and Z — > pp analyses, correction factors to account for the efficiency and detector resolution 

166 within the selected acceptance based on the simulation are calculated differentially in event centrality, p 2 , and y 2 . 

167 In each decay channel, the correction factor is applied and the background, estimated by the same-sign charged 

168 pairs, is subtracted. The two decay channel measurements are then combined with weights set by their respective 

169 uncertainties. The fully corrected p 2 distribution for the Z boson is shown in Fig. [3] The data distributions plotted in 

170 (0-80)% centrality agree in shape with the Pythia simulations of Z boson production in p + p collisions. 

171 The nuclear modification factor, R^c, normalized to the statistically-robust central collisions bin can be defined as: 

R = <AU)(C) (l/W evt , P )d 2 AWd3;dp T 
PC (N con )(P)(l/N^c)d 2 N c /dydp T 
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Figure 3: Left: Efficiency corrected yields of prompt photons in \rj\ < 1 .3 in the isolation cone radius R lso = 0.3 and isolation energy of 6 GeV, 
scaled by (Taa) from Ref. 1141 . Statistical errors are shown by the error bars. Systematic uncertainties on the photon yields are combined and 
shown by the yellow bands. The scale uncertainties due only to (7*aa). The solid lines are Jetphox calculation of photon yields and dashed line 
is the Pythia. Corrected (middle) and rapidity (right) distributions of the measured Z — » l + l~ (I = e,/i) data compared to Pythia predictions 
normalised by area from Ref. 1111 . The error bars represent statistical uncertainties, and the filled bands represent systematic uncertainties. The 
lower panes display the ratio of data over Pythia. 



172 where (N co n)(F) and (A^ co u)(C) are the number of binary nucleon-nucleon collisions calculated for peripheral and 

173 central collisions respectively, and (l/N evt c)d 2 Nc/dydpT and (1 /N ew tp)d 2 Np/dydpT are the differential yields of Z 

174 per-event for central and peripheral collisions respectively. 

175 The nuclear modification factor for the Z boson as a function of pj is shown in the left panel of Fig. Hfor several 
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Figure 4: Left: the Rpc of Z boson for four centrality bins. The error bars show the statistical uncertainties, and the shaded areas at the point 
location, the residual systematic uncertainty from Ref. 1 1 1 1. Additional uncertainty, coming from the (A' co ii)(P) to (N co \i)(C) ratio is the same for 
all points in the panel and is shown centered about unity. Right: photon yields in \r]\ < 1.3 in isolation cone radius i? lso = 0.3 and isolation energy 
of 6 GeV, divided by Jetphox 1.3 predictions from Ref. [14]. Statistical errors are shown by the error bars. Systematic uncertainties on the photon 
yields are combined and shown by the yellow bands. 
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176 centrality bins. The systematic uncertainties largely cancel out in the ratio, however the statistical uncertainties remain 

177 large, especially for the most peripheral (40-80)% bin. The statistical uncertainties are correlated in different panels 

178 of the figure, because all ratios use the same (0-5)% centrality data as the denominator. The right panel of Fig. [4] 

179 shows the analogue results for isolated photons presented as ratios of the photon spectra divided by the results of the 
iao predictions from the Jetphox and by the (7aa)- Within the uncertainties, the Rpc of Z bosons, and the ratio of photon 
iai production to the prediction values are consistent with unity. 
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Figure 5: W boson Rpc as a function of collision centrality from Ref. |5 |. The error bars are uncertainties from the fit method and statistical 
uncertainties. Yellow bands include uncertainties from the number of binary collisions. The dashed line is the result of a flat line fit. 



182 Figure [5] shows the Rpc of the W boson as a function of centrality. Within the uncertainties, which are mainly 

183 statistical it remains consistent with unity. 

The left panel of Fig. [6] shows the yields of isolated photon scaled by (Taa) as a function of (N p - drt ). The Z boson 
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Figure 6: Left: centrality dependence of Z boson yields divided by (N co \{), measured in \y z \ < 2.5. Results for ee (upward pointing triangles) and 
/j/i (downward pointing triangles) channels are shifted left and right respectively (for visibility) from their weighted average (diamonds) which 
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for each selected centrality interval. Statistical errors are shown by the error bars. Systematic uncertainties on the photon yields are combined and 
shown by the yellow bands. 
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184 yields per-event, divided by (N co \\), are shown in the right panel of Fig.[6]for Z — > ee and Z — > //// independently. The 

185 figure demonstrates that the two channels are consistent within their uncertainties and can be combined together. The 

186 combined results are also shown in the right panel. All values are integrated over the pseudorapidity range in which 

187 the corresponding particles are measured. The results for photons and Z bosons are presented in several p T slices 

188 and for W and Z bosons integrated over all pj. Figure [6] shows that the yields of all measured bosons scale with the 

189 number of binary collisions {N co u) (or nuclear thickness function (7aa)) calculated with the Glauber model. 

190 4. Conclusions 

191 The ATLAS experiment measured the W, Z boson and isolated photon production in Pb+Pb collisions at y%JV= 

192 2.76 TeV using data collected in the 2010 and 2011 LHC physics runs. Isolated photons are reconstructed in the 

193 range p T = 45 - 200 GeV and for \t]\ < 1.3. The W boson yields are measured via single muons registered in 

194 the range \rj\ < 2.5 and with pj > 7 GeV. Z bosons yields are fully corrected to \y z \ < 2.5 for the mass region 

195 66 GeV < mz < 116 GeV. Analyzed Z — > ee and Z — > pp decay modes produce consistent results. The measured 

196 distributions of isolated photons agree in shape and in yield with the Jetphox predictions calculated without taking 

197 into account nuclear modification. The shape of the pj and y distributions for Z are consistent with the prediction of 

198 Pythia generator. The total yields of all measured bosons scale proportional to the nuclear thickness function, (Taa) 

199 (or {N co \{}). The Nuclear Modification Factor Rpc for Z and generator normalized for photons agree with unity as a 

200 function of pj. 
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